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We show that universal quantum control of a two-electron singlet-triplet spin qubit can be achieved 
using Landau-Zener-Stiickelberg interferometry. Going beyond normal Landau-Zener dynamics with 
infinitely long constant velocity sweeps across an energy level anti-crossing, we focus on a physical 
system consisting of a two-electron double quantum dot, where the spin states can be admixtures of 
charge states and the level velocity can be tuned in a time-dependent fashion. Our results indicate 
that charge coherence must be treated on an equal footing with spin coherence. In particular, 
we predict the presence of finite-time effects, which result in population transfer even in cases of 
an incomplete sweep through the anti-crossing. The competing requirements of adiabaticity and 
coherence are reconciled using specially designed pulses with a tunable level velocity. As a relevant 
example, we demonstrate that a Ifadamard gate can be implemented for a realistic set of conditions 
in a GaAs double quantum dot device. 



I. INTRODUCTION 

Electron spins trapped in quantum dots (QDs) are 
promising candidates for implementing a scalable quan- 
tum computer^ ^. Most of the DiVincenzo criteria have 
been achieved with spin qubitSjincluding initialization'^, 
readoulP', and coherent controF^. While experimental 
progress has been impressive, many of these methods are 
not yet accurate enough to allow large scale quantum 
computing with single spins. A key challenge can be 
appreciated by considering the relevant timescales asso- 
ciated with the spin dynamics. In GaAs quantum dot 
devices, it is well known that the hyperfine interaction 
leads to a randomly fluctuating nuclear field, i3n /"^ 2 
mT, which results in a 10-20 ns spin dephasing tim d^l^l^l 
For comparison, the Rabi period obtained in a GaAs dou- 
ble quantum dot (DQD) using conventional electron spin 
resonance (ESR) is on the order of 110 n^. The use of 
spin-orbit driven electric dipole spin resonance (EDSR) 
in GaAs leads to even slower Rabi periods of roughly 210 
nfl. Therefore, in the case of single spin rotations, gate 
operation times are nearly an order of magnitude slower 
than the inhomogeneous spin dephasing time. This issue 
is not specific to GaAs based nanostructures. In InAs 
nanowires, where the spin-orbit interaction is larger than 
in GaAs, a Rabi period of ^ 17 ns was reported. How- 
ever, it still is twice as long as the spin dephasing time 
of ~ 8npl. 

Viewed from a different perspective, the maximum ac 
field generated in DQD ESR experiments is on the order 
of 2 mT, which is the same magnitude as the fluctuating 
nuclear fielcP. As a result, single spin rotations in GaAs 
qubits follow imperfect trajectories on the Bloch sphere, 
resulting in a reduced oscillation visibility and gate er- 
rord^. Single spin selectivity imposes an additional chal- 
lenge to the development of a spin-based quantum pro- 
cessor. Magnetic fields generated in ESR are difficult to 
localize on the nanometer scale. Without g-factor con- 
trol, or local magnetic field gradients, the spins located 



in a quantum register would rotate at the same rate in 
the presence of a global ESR field. The long term goal is 
to be able to drive selective single spin rotations, with- 
out affecting neighboring spins that are on average only 
20-50 nm away. 

Instead of using the spin-up and spin-down states of a 
single electron, the qubit basis states can be represented 
by two (out of four) two-electron spin states confined in 
a double quantum dot (DQD)^^. For a qubit whose basis 
states are encoded in the singlet S and triplet Tg spin 
states of a DQD, the two-electron exchange interaction 
allow s for a fast 200 ps VSWAP gat(P. Recent experi- 
ment Ji^'^ have also demonstrated the possibility of real- 
izing a conditional two-qubit gate. Within this two-spin 
version of a qubit, the two-qubit gate realized in Ref. [51 
can be interpreted as a single qubit operation. However 
the exchange interaction only allows for rotations about 
a single axis, whereas to generate arbitrary rotations one 
needs two perpendicular rotation axes. The generation 
of a nuclear magnetic field gradientP^ provides rotations 
about a second, non-collinear axis. While remarkable, 
this method also presents some difficulties when it has 
to be extended to a large number of qubits. It requires 
that the nuclear polarization be controlled in each DQD 
to create the desired gradient field. An advantage of this 
method is that the generation of the nuclear field gra- 
dient reduces nuclear spin fluctuations, resulting in an 
increase in the spin dephasing time^^. 

Recently, it has been proposed to use a two-spin basis 
consisting of the singlet S and triplet T+ spin stateJ^SHUl, 
Quantum control of the S-TV qubit relies on Landau- 
Zener-Stiickelberg-Majoran£p2Hlll (LZSM) physics, which 
occurs in the system when the S-T-|_ qubit is repeat- 
edly swept through the hyperfine mediated S-T+ anti- 
crossing. This method, in addition to being all-electrical, 
has the advantage of being local. 

LZSM physics describing the passage of a two-level 
quantum system through an anti-crossing can be ap- 
plied to different fields of physics and chemistrj{2ll. In 
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quantum information science, LZSM theory accurately 
describes the observed interference fringes (Stiickelberg 
oscillations) of a superconducting qubit driven back and 
forth through its anti-crossing^, as well as the coherent 
manipulation of a two-spin qubit encoded in the S and 
triplet T+ spin states, with dynamics driven by rep eated 
passages through a hyperfine-mediated anti-crossin ^^^ l ^^ l 
In this article we develop a master equation to de- 
scribe partial adiabatic passages between the spin singlet 
S and triplet T_|_ states at the hyperfine mediated anti- 
crossing in a S-T+ spin qubit. With our new formalism 
we show that charge dynamics can significantly hinder 
LZSM interferometry of spin states. While most of the 
interesting spin dynamics happens in the (1,1) charge 
configuration, with one electron per dot, initialization 
and measurements are done deep in the (2, 0) charge con- 
figuration, where both electrons are in the left dot. Here 
{I, r) denotes the number of electrons in the left and right 
dot. Crossing the (1,1) O (2,0) interdot charge transi- 
tion necessarily involves charge dynamics. Since charge 
states dephase on shorter time scales than spin states, 
it is essential to consider this fast effective de cohere nce 
mechanism while studying spin qubit evolutiorPSl2l!. In 
particular, we use our formalism to analyze spin-charge 
dynamics associated with detuning pulses that have a 
tunable level velocity, with a high level velocity away 
from the anti-crossing and a slow level velocity in the 
vicinity of the anti-crossing^^. Furthermore, we show 
that this type of pulse is suitable to realize a Hadamard 
operation, given by the single-qubit unitary 
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with a gating time smaller than the characteristic inho- 
mogeneous dephasing time Tj. The Hadamard gate is 
a key element of many quantum algorithms. It is used, 
for example, in Shor's quantum algorithm for factorizing 
large integer^^H, 

The paper is organized as follows. We start by review- 
ing standard Landau-Zener transition theory, which is 
valid for an infinitely long ramp through an anti-crossing 
with a constant level velocity. In Sec. [H] we review the 
solution of the finite-time LZSM model, which can be 
used to model realistic experiments, and we demonstrate 
within the scope of this theory how fine tuning of the 
level velocity can be used to increase the visibility of the 
quantum oscillations. Section HI focuses on the phys- 



ical implementation of LZSM physics in a two-electron 
spin qubit. We derive an effective Hamiltonian describ- 
ing the dynamics of the states in the vicinity of the S-T+ 
anti-crossing. Compared to previously derived effective 
HamiltonianJl^, we include the effects of a charge admix- 
ture in the spin-dependent anti-crossing. The last part of 
the section is devoted to the derivation of a master equa- 
tion that describes the evolution of the density matrix. 
In Sec. |IV| we show the results for the singlet return prob- 
ability obtained from the solution of the master equation 
and explore the effects of a time-dependent level velocity 



and charge induced decoherence. 



II. ADIABATIC CONTROL OF A QUANTUM 
TWO-LEVEL SYSTEM 



There are numerous problems in physics that deal with 
the physics of two-level systems. The most common ex- 
ample is Rabi's formulcpS which describes the occupa- 
tion of a two-level system that is excited by a coher- 
ent field. In quantum information science, Rabi oscil- 
lations are widely used to manipulate an electron spin 
confined in a QEPH, Another widely studied problem 
involving only two quantum levels is adiabatic passage, 
which is commonly employed in nuclear magnetic reso- 
nanc^^. The physics of adiabatic passage can be found 
in a variety systems, and several theoretical models have 
been developed to describe different kinds of adiabatic 
processepilt^. There is however a particular model that 
has proven to be applica ble in many distinct fields of 
physics, the Landau-Zenei'^^'^ni model. We refer to it in 
this article as the LZSM model, since it was indepen- 
dently studied by Stiickelber and MajorangP^in 1932. 

The Hamiltonian studied in the LZSM model describes 
a system with two energy levels [see Fig. [IJa)] that are 
coupled by an off-diagonal matrix element X, H = —eaz 
+ Xcr^. Here e is a detuning parameter that is controlled 
experimentally. The LZSM model considers a detun- 
ing sweep described by a constant level velocity, a — 
d(i?2 — Ei)/dt, i.e., e = at/2. The main result of the 
theory is the asymptotic expression for the non-adiabatic 
transition probability when the propagation lasts from 

ti — —OO to t{ — CXI, 
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Here we use a generalization of the LZSM model, 
known as the finite-time LZSM modeP^I. It resolves the 
problem of the energy divergence when j = ±oo, and in 
contrast to the simple case it gives the relative phase be- 
tween the states, which is a crucial for predicting the time 
evolution of any quantum system. Particularly, knowl- 
edge o f the relative phase is essential in LZSM interfer- 
ometrjEES! jn which the system is driven back and forth 
across an anti-crossing. The driving generates quantum 
interference between states which are directly observable 
in the non-adiabatic (or adiabatic) transition probability. 



A. Finite-time 
Landau-Zener-Stiickelberg-Majorana propagator 

The unitary evolution operator defined by the LZSM 
HamihoniaiPSMl^ 
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Figure 1. (color online) (a) The Landau-Zener-Stiickelberg- 
Majorana (LZSM) problem. The LZSM model addresses the 
problem of a two-level system that is swept through an anti- 
crossing. It assumes an infinitely long ramp (i.e. propagation 
from ti — >■ — cxD to tf — >■ oo), a constant coupling constant A 
leading to a splitting 2 j Aj at t = 0, and an energy difference 
between the levels which is linear in time AE{t) = at. Its 
main result yields the non-adiabatic transition probability, 
ftzsM. An extension of the model to finite-times resolves 
the problem of infinite energies and undefined phases, (b) 
Singlet and triplet T+ energies in a DQD as a function of 
detuning e, where ei — e{ti) and £f — e{tf). The two-electron 
DQD is a physical realization of the LZSM model. Here, the 
hyperfine coupling of the two spin states leads to a splitting 
Ahf. (c) The LZSM model assumes a constant level velocity 
a. In order to increase adiabaticity one can lower alpha, but 
in order to keep short pulses it is preferable to use multi rise- 
time pulses, (d) Comparison of possible pulses that can be 
used to manipulate the singlet triplet T+ qubit. In terms of 
total propagation time, "double hat" pulses represent a good 
compromise between conventional trapezoid and convoluted 
pulses. 
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+ Di^2 (e'^Tf^j Di^2 (^e * r;^ 



Here t = ^ajht is a dimensionless time, r/ = X/\/ah is 
a dimensionless coupling, T{z) is the gamma functioiP^, 
and D^{z) is the parabolic cylinder functioiP^. 

The LZSM propagator fully determines the partially 
adiabatic dynamics of a quantum two-level system. In 
the case where the two-level system encodes a qubit, 
Eq. ([4]) allows the design of single qubit operations. 
This method has been used to control superconducting 
qubits^"^ and more recently to imple ment a qubit en- 
coded in the spin of a two-electron stat c!^^ * ^^ -*. However, 
since the spin states are weakly coupled, < 1, it is 
hard to achieve an equally weighted coherent superpo- 
sition of spin states and fully explore the entire qubit 
state space. In order to achieve full control over the spin 
qubit, it would be necessary to increase "q by making a 
smaller. However, the LZSM equation requires an ex- 
ponential increase in the propagation time in order to 
achieve a fully adiabatic transition. In a real physical 
system, this method is unpractical because the pulse du- 
ration should remain shorter than the coherence time of 
the two-level system. 



B. Enhancing Adiabatic Control 

As demonstrated in Ref. [TSl it is possible to use more 
complex pulses to increase the balance of the populations 
while keeping the manipulation time below the decoher- 
ence times. The key idea relies on an observation based 
on the finite-time LZSM model. For a slow level velocity 
a, which favors adiabatic passage, most of the popula- 
tion change occurs in the vicinity of the anti-crossing. It 
is therefore possible to use detuning pulses that have a 
time-dependent level velocity. Let us consider two types 
of pulses, as illustrated in Figjljc) and (d)). The first is 
a conventional linear pulse, which is standard in LZSM 
theory. The second pulse profile consists of linear detun- 
ing ramps in a fast-slow-fast rise-time sequence, which 
we refer to as multi-ramp pulse. The unitary evolution 
of such a general sequence can be written using Eq. @ 
as 



U{t{,ti) = (7fast2(if, t2)Uslowit2, il)C^fastl (^1 , ^i) 

ti) Ui2itl,ti)\ (7) 
ti) U22{ti, ti)) ' 
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Here ti = thi, t2 = t^i + tsr, and tf = t^i + tsr + ifr2, 
where tf^j is the propagation time associated to the jth 
fast sequence, and tgr corresponds to the slow sequence. 
We use this notation to refer to the corresponding level 
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Although the passage through the anti-crossing is done 
under conditions that largely favor non-adiabatic tran- 
sition, the multi-ramp pulse leads to a more adiabatic 
passage than a linear pulse with the same total propaga- 
tion time At — tf — ti. 

All the interest of the method lies on this last property, 
if we think of a two-level quantum system as a phys- 
ical implementation of a qubit. Since a qubit cannot 
be perfectly isolated from its environment, any coherent 
state is going to decohere within a characteristic time 
scale. In the case of LZSM driven qubits, this implies 
that there must be a compromise between the optimal 
time to realize a quantum operation and the decoherence 
time. However, the use of multi-ramp pulses allows an 
increase in adiabaticity while keeping propagation times 
smaller than decoherence times. 

This control method is more efficient in the case of 
strong coupling, 77sr > 1- It allows to achieve fast full 
population transfer by fine tuning of ti (or t[). We present 
in Fig. [2jb) results obtained for Pna as a function of ti 
for A = 5 /ieV and t2 — ti — 4.0 ns. A method allowing 
for high fidelity full population transfer using composite 
chirped pulses has also been proposecP^. 



Figure 2. (color online) (a) Comparison of the non-adiabatic 
transition probability, Pna, for a "double ramp" (orange) and 
a linear pulse (red) as function of ti. Here, t\ — ti = 0.1ns, 
t2 — ti = 2.0 ns, if — t2 = 0.1ns, A = 50neV, Qfri — 3 • 
lO'^ eVs-\ a,r = 5 • lO"* eVs-\ and om = 2-10'^ eVs'^ The 
condition that the system must reach the same difference in 
energy defines the level velocity for the Gaussian pulse. The 
double hat pulse sequence increases the adiabatic transition 
probability, while keeping the total propagation time fixed, 
(b) In the strong coupling regime the double hat pulse allows 
full adiabatic population transfer by fine tuning ti. Here A = 
5 neV and t2 — ti = 4.0 ns, the others parameters are the same 
as in (a), (inset) Schematic representation of the pulses used 
to obtain non-adiabatic probabilities. The colors of the pulses 
are matched to the results shown in (a) and (b). 



III. ADIABATIC CONTROL OF A S T+ QUBIT 

In the following, we apply the previously developed 
ideas to a physical implementation of a LZSM driven 
qubit. We focu s on t he two-spin S — T+ implementation 
in a GaAs DQEffSEl^ see Fig.|l];b). Although the dynam- 
ics of the system under multi rise-time pulses has already 
been studied experimentalljff^, there is still a need to gain 
a better understanding of the charge-noise-induced spin 
dephasing. We will show, among other things, that the 
measurement of finite-time LZSM oscillations can pro- 
vide a tool to qualitatively define the strength of charge 
noise. 



velocities aj, dimensionless times Tj, and dimensionless 
couplings r/j. 

To illustrate the advantage of a multi-ramp pulse, we 
plot in Fig. [2] the non-adiabatic transition probability, 
^na = 11*11(^1, tf) I , as a function of ti for two differ- 
ent coupling constants. Here a non-adiabatic transition 
refers to a transition where the system remains in the ini- 
tial uncoupled energy eigenstate. The level velocity for 
the linear pulse, agi is obtained by the condition that the 
system must be detuned to the same difference in energy, 

afilifrl + asr^sr + Q!fr2tfr2 /o^ 

"si = • (8) 

tf — ti 

We start considering a weak coupling regime for the slow 
component of the pulse, r/sr < 1- The results are pre- 
sented in Fig. [2ja), where we have chosen ti —ti = 0.1ns, 
t2 — ti = 2.0 ns, tf — t2 = 0.1ns, A = 50neV, af^i = 
3-lO^eVs-i, = S-lO^eVs-^, and af^ = 2-106eVs-i. 



A. Double Quantum Dot Spin States 

The spin preserving part of the Hamiltonian describing 
the confinement of electrons in a DQD in the presence of 
an external magnetic field can be written using a simple 
two-site hopping model, 

4=1.2 ^ ^ i 

+ (cLc2a +h-c.) . 

(9) 

The index j = 1, 2 labels the dot number and a =t,|= 
±1 the spin of the electron. We denote the energy of 
a single confined electron by and the Zeeman energy 
associated with its spin is given by g* fiBBa /2, where g* 
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Figure 3. (color online) (a) Energy diagram for the relevant 
states in the DQD as a function of e. The spin states for the 
implementation of the qubit are the hybridized singlet S and 
the triplet T+. 



denotes the effective Lande g-factor, ji-g. the Bohr magne- 
ton, B the strength of the external magnetic field. The 
operators Ci^ and describe respectively the annihi- 
lation and creation of an electron in dot i with spin a. 
Two electrons occupying the same QD cost an intra-dot 
Coulomb energy u. The last term of Eq. ^ accounts 
for electron tunneling between the dots with strength r. 
We neglect the inter-dot Coulomb interaction because it 
only produces a constant shift of the energy levels. 

Since most recent experiments on a DQD system are 
operated in a regime with at most two electrons, we can 
project Eq. (|9| into the subspace sp anned by th e charge 
configurations (0,2), (2,0), and [i^i^M^m, The di- 
agonalization of the resulting Hamiltonian leads to six 
low-energy states which are superpositions of the sin- 
glets S(0,2), S(2,0), and S(l, 1) as well as the triplets 
To(l, 1), T+(l, 1), and T_(l, 1). The triplet states with 
two particles in the same dot must have electrons occupy- 
ing higher-energy orbitals levels due to the Pauli princi- 
ple. This results in a relatively high separation in energy 
(^ 400 /icV) from a singlet with both electrons occupy- 
ing the same dolP^. Consequently the triplets states with 
two electrons in the same dot can be safely neglected for 
the purpose of the current study. 

An energy level diagram is shown in Fig. [3j where we 
plot the energies of the relevant two-electron spin states 
as a function of detuning e = ei — £2- The degeneracy 
of the singlets S(l, 1) and S(2, 0) at e = u, as well as 
between S(l,l) and S(0, 2) at e = — u, is hfted due to 
tunneling, which results in a splitting of the levels by 
2^/2T. The de generacy between the spin singlet S and 
spin triplet Tq is lifted due to the exchange interaction^®. 
This property has allowed the encoding of a spin qubit 
in these two spin states and its manipulation it via the 
exchange interaction^ . 

Here we concentrate on a particular value of detuning 
that we have denoted by Ec in the energy level diagram. 
This point corresponds to the crossing of the singlet state 
with the triplet T+ state. What renders this point special 
is that in reality there is an anti-crossing due to the hy- 



perfine interaction between the electron spins and nuclear 
spins of the host material. It has been demonstrated ex- 
perimentalljff^ and theoreticalljEl that coherent control 
of the singlet S and triplet T+ can be achieved by detun- 
ing the system from an initially prepared S(2, 0) through 
the hyperfine mediated anti-crossing. 

The hyperfine interaction is described by the Hamilto- 



nian 



Hyiv — Si ■ hi + S2 ■ h2 



(10) 



between the electron spin Si and the effective magnetic 
fields hi that are generated by the nuclear spins li in 



dot 



k rk 



Afl. 



I. The Overhauser field operators hi = X^fcLi 
describe the nuclear bath. Here rii is the number of nu- 
clei in dot i and — Wifcfo |'!/'i(''fc)l is the hyperfine 
coupling constant with the fc-th nucleus in dot i, with 
''pii'fk) the electron wave function, vq the volume of the 
unit cell and Vik the hyperfine coupling strength. A more 
convenient form of the hyperfine interaction is obtained 
by introducing the spin ladder operators Sf" = Sf ± \Sf 
and = hf ± i/if , which yields 



HF 



2^ 



{2S!ht + stK + s-ht). (11) 



The longitudinal part of Hup is diagonal and its contri- 
bution will add to the energy of the triplet state. The 
transverse part 



^HF 



(12) 



generates the so-called fiip-flop process that result in an 
energy gap at Ec and allow for mixing of the S and T+ 
spin states®^ . 



B. S — T+ Effective Hamiltonian 

In order to derive the Hamiltonian describing the dy- 
namics of the two-electron spin states near the hyperfine 
induced anti-crossing, we consider a 3 x 3 Hamiltonian 




(13) 



which accounts for the dynamics of the T+(l, 1), S(l, 1), 
and S(2,0) spin states in an external magnetic field, B. 



The diagonalization of Eq. ( 13 ) yields two hybridized 



singlets and the triplet state T 



|S)=c(£)|S(l,l)) + Vl-c(e)2|S(2,0)), 

IS') =c'(e)|S(l,l)) + v/l-c'(£)2|S(2,0)), 
|T) = |T+(1,1)), 



(14) 

(15) 
(16) 
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with respective energies 



The charge admixture coefficients c(e) and c'(e) are 
e — u — \/St^ + {u — e)2 



cie) = 
c'ie) = 



^8t^ + i-e + U+ V8r2 + (£"3^)2 
e-u + a/8t2 + - £)2 



(17) 
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We can now make the following basis change 
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T(e) - c{e) c'je) 



(21) 
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such that 



Figure 4. (color online) The effective hyperfine mediated cou- 
pling between |S) and |T) is a function of detuning. This 
is a consequence of the state |S) being an superposition of 
different charge states. 



The longitudinal part of -ffuF can be included in the 
energy of the triplet state g* hbB g* ^b{B + B^^ + 
B^2): while the transverse part 



Ht^{e)^T{s)^H,{s)T{e) 

^g*fiBB 
£;s 
Es'. 
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^1 Btt)- 



(27) 
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mixes the spin states |S) and |T). We therefore find that 
the full Hamihonian, in the {T+(l, 1), S(l, 1), S(2, 0)} ba- 
sis, describing the dynamics of the singlet S and triplet 



We add now to the Hamiltonian defined in Eq. (flSl) the T+ near the hyperfine anti-crossing is given by 



hyperfine interaction defined in Eq. (10 1. Since we are 
not interested in describing the nuclear spins dynamics, 
but its effect on the two-spin states S and T+, we can 
model the action of the Overhauser field operators by 
introducing a classical stochastic variable whic h acco unts 
for fluctuations in the nuclear spin ensembl d^ * ^ * ^^ l By 
setting 



ffs-T+(e) 




(28) 



where B = B + Sf^ ^ + S,^ 2 and 



hi = g*fJ.B'Bn,; 



(24) 



we can interpret ^n,i as the effective random magnetic 
field interacting on Si. Under normal experimental con- 
ditions we have k^T 3> (?nMn^, where g^ and /in are the 
nuclear g-factor and magneton. The nuclei can, in this 
limit, be assumed to be completely unpola rized, resulting 
in a Gaussian distribution of nuclear field^^'^ESl 



A = (S(l, 1)\hUt+{1, 1)) = g*fiB{B-^ - B-,)/2V2. 

(29) 

In the basis that diagonalizes Eq. ([I3|, the Hamiltonian 
defined in Eq. ( 28 1 reads 



ffs-T+(e) 



'g*piBB c(e)A 
c(e)A Es(e) 
d{e)\ 




(30) 



P{Bn,t 



1 
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(25) 



with a = A/ g* ^^y/n and the hyperfine coupling constant 
v4 « 90 /ieV, and the approximate number of nuclei over- 
lapping with the electronic wave function n « 10^ — 10^. 
By defining B^ ^ = i?^ ^ ± j, the hyperfine Hamilto- 



nian can be written by analogy with Eq. ( 11 ) as 



The projection of Eq. (|30|) onto the Hilbert space 
spanned by {|S), |T)} leads to the effective 2x2 Hamil- 
tonian describing the dynamics at the singlet-triplet T+ 
anti-crossing. The higher energy singlet state |S') can 
be safely neglected if the pulses bringing the states |S) 
and |T) into resonance are fully adiabatic with respect 
to the energy separation AiJgg/ — 2\/2t of |S) and |S') 
at the charge anti-crossing. Taking into account that the 
detuning is time-dependent, e — s(t), we finally obtain 



^HF = J E C^^^K^ + S^B-, + S-Bl^) . (26) H{t) 



Es{e{tm{S\+Err\T){T\+g{e{t)) (|S)(T| +h.c) , 

(31) 



7 



where Eg is defined in Eq. (17 1 and Et = g*iiBB. 
This effec tive Hamiltonian differs form previous deriva- 
tion d^^ l ^^ in the coupling g{t) — c{e{t))X, which is time- 
dependent. Here, c(e(t)) and A are given in Eqs. (20) 
and (29). As it can be seen from the functional form 
of (7(t)7the effective coupling strength between the spin 
states depends on the charge state (Fig. [4]). This re- 
sult is rather natural since the matrix element between 
S and T+ goes to zero when the detuning is such that 



(S|7?|Lp|T+) 







S = S(2,0) and T+ = T+(l,l), 
for e ^ u. It also implies that the physics described 
by Eq. ( [sT] ) goes beyond standard LZSM theory with a 
constant coupling. However, since both Hamiltonians de- 
scribe adiabatic passage through an anti-crossing, we can 
assume that the dynamics are qualitatively similar. We 
can therefore expect that our previous discussion about 
enhancement of adiabaticity based on LZSM physics re- 
mains valid. Furthermore, due to the peculiar form of 
g{t) (c.f. Fig. |4]), we can assume that it is possible to ob- 
serve finite-time interferometry phenomena in close vicin- 
ity of the anti-crossing. 



C. Master Equation 

In order to compare our theory with experimental mea- 
surements, it is not sufficient to solve the Hamiltonian 



dynamics provided by Eq. (31 ), because some important 



phenomena that can influence the outcome of the exper- 
iment are not taken into account. Among these is the 
spin-orbit interaction that leads to phonon mediated re- 



laxation and charge fluctuations that lead to dephasing 
of a charge superposition state. These phenomena can 
be taken into account in a master equation formalism for 
the density matrix. In the Lindblad formalism, the time 
evolution can be expressed as 



Ljp,L] 



L„pL^i\)- (32) 



The operators Lj are called Lindblad operator^^^EI] and 
they describe the dissipative effect of the environment 
on the system in the Born-Markov approximation, which 
consists of two assumptions that lead to Eq. (32). The 



Born approximation supposes a weak coupling between 
the system and the bath, in the case of spin-orbit inter- 
action to the phonon bath, while the Markov approxima- 
tion consists in neglecting any type of memory effects of 
the bath during the system evolution. To fully describe 
the effect of the environment, one needs — 1 opera- 
tors where N is the dimension of the system's Hilbert 
spacd^215ij -ppj. ^]-^g j,g^gg q£ ^ two-level system, the Lind- 
blad operators are Li — -^/tTct-, L2 = •\/r+cr+, and 
L3 — yJT^pOz^ where (t_ and oj^ are spin ladder oper- 
ators and (T2 is the z Pauli matrix. They respectively 
describe relaxation from the excited state to the ground 
state with rate r_, relaxation from the ground state to 
the excited state with rate r+, and pure dephasing with 
rate F^. 

By substituting Eq. 



( |3l| into Eq. ( 32 ) and using the 
expression of the Lindblad operators, the first order dif- 
ferential equation for the S-T-|_ density matrix can be 
written as 





-^^ 


P12 _ 




P21 1 


-bit) 


\p11) 







-bit) 



-4F^) 



-bit) 



^t) - \ {T, 

bit) 



r_ + 4F„ 



-bit) 
bit) 
-r- / 



\ 






P12 








\P22/ 



(33) 



where we denote the uncoupled ground state and excited 
states with the indices 1 and 2. This system of four 
coupled ordinary differential equations can be reduced to 
a system of three coupled ordinary differential equations 
(Bloch equations) by introducing new variables defined 

by 



X = P12+ P21, 

y = i(Pl2 - P21) : 
Z = Pll- P22- 



(34) 



This set of variables is completed by the conservation of 
probability condition. 



Pii + P22 = 1- 



(35) 



Substituting Eqs. (34) and (35) into Eq. (33), one finds 



the system of ordinary differential equations for the new 
variables, 



A-- ^' ^^^ y-^iT++r^)x-2r^x, (36) 

Es - Et- a 1 
y= \ ^ a:-2|z--(F++r_)t/-2F^y, (37) 



i = 2|y-(F++F_)z + F_-^^ 



(38) 



Since F^ and F_ can be related to each other by consid- 



ering the limiting case of thermal equilibrium, Eqs. (36) 



(38) can be simplified to include only two independent 



rates. If the system reaches thermal equilibrium then the 
detailed balance equation pflT^ — P22^~^ holds. More- 
over, the populations are given by the canonical ensem- 
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Figure 5. (color online) (a) Relaxation rate ri(£) obtained for 
B — 100 mT. Spin-orbit mediated flip-flops of the electron 
spin accompanied by emission or absorption of phonons is 
assumed to be the leading relaxation mechanism. The relax- 
ation rate is maximal for £ = Sc where |-Bs(s) ~ ^ ksT. 
This results in a strong mixing of the spin states due to ther- 
mal fluctuations, (b) Charge noise induced dephasing rate 
r2(£). The dephasing rate is assumed to have a functional 
dependence proportional to 1 — c{e)^ to ensure that dephas- 
ing related to a charge noise is suppressed when the system 
is in a (1, 1) charge configuration. 



ble, plf — exp{—PEi)/Z, where /3 — l/k^T, where /cr is 
Boltzman's constant, T is the temperature, and Z is the 
partition function. This yields 



Pll 



(39) 



where we used hio 



TT'd I 
T 



The energies E§ and 



are obtained by diagonalizing Eq. (31 1. 
Since F- is associated with a process involving the 
emission of a phonon, we can express it as 



F_ =7i(l + n(w)) 



(40) 



where n{uj) is the Bose-Einstein occupation number, 
n{uj) = {exp{/3fiuj) — 1)^^, and 71 is the spontaneous 
decay rate, which we treat as a phenomenological pa- 
rameter here. Similarly we obtain 



(41) 



Combining the results of Eqs. (39 1, (40), and (41), we 
obtain the relaxation rate 



ri(e) = 



71 coth 



\Ei{s)-E^{s) 



(42) 



This function is plotted against e in Fig. [5^, and has a 
peak at £ = where \E^{e) — E^ \ <^ k^T resulting in a 
strong mixing of the states due to thermal fluctuations. 

Pure dephasing can originate from orbital effects. In 
our current description of the problem, we have ne- 
glected that the wave functions of the singlet |S(2, 0)) and 
triplet |T+(1,1)) couple differently to the background 
charge environment due to their different charge con- 
figurations. Thus, a superposition state of the form 
= q;|S(2, 0)) -I- /3|T+(1, 1)) is sensitive to background 
charge fluctuation s (cha rge noise), which leads to dephas- 
ing of the state If the qubit is in a superposition 
with the same charge state, \lp) = a|S(l, 1))-|-/3|T+(1, 1)), 
then the effects of charge noise are negligible. Therefore, 
the dephasing can be written as 



r^(e) = 72(i-|c(£)r). 



(43) 



where 72 is the charge noise rate. This form ensures that 
charge noise induced dephasing is suppressed when the 
spin states have the same charge configuration. In Fig. [5] 
we plot the rates ri(e) and T^{e). 

Finally, the Bloch equations, written in a matrix form, 
and describing the dynamics around the S-T+ anti- 
crossing are 








4r^(£W)) 



~Es{e{t))+ET 
\h{Ti[e{t)) + 2T^{e{t))) -2\c{e{i)) 



2\c{e{t)) 



-hT^{e{t)) 




±h 



(44) 



The different signs in front of the inhomogeneous term 
come from the fact that the states |S) and |T) exchange 
their roles as ground and excited state of the system at 
e = £c- 



IV. RESULTS 

We now compute the singlet return probability for 
two different detuning functions e{t). We first consider 
the case where £{t) is obtained through convolution of 
a trapezoid pulse with a Gaussian. We then consider 
what we call a "double hat" pulse [see Fig. [l]jd)]. This is 
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a pulse consisting of a fast-slow-fast rise-time succeeded 
by a waiting time and a reverse fast-slow-fast rise-time. 
Both types of detuning pulses have already been experi- 
mentally used to demonstrate control over a S - T+ based 
qubitie'is. 

As stated previously, here we are interested in studying 
the effect of charge dynamics on the partial adiabatic 
spin dynamics in the vicinity of the hyperfine mediated 
anti-crossing. We therefore solve Eq. (44) for different 
values of 71 and 72. The singlet return probability is 



tw (ns) 



obtained by solving Eq. ( 44 ) for different realizations of A 



that obey the probability distribution defined in Eq. ( 25 1 



Since the actual number of nuclear spins interacting with 
the electron spin is unknown, we use the experimentally 
measured hyperfine splitting from Ref. [16] to determine 
the standard deviation. The average value of Ps is then 
calculated according to 



1 ^ 1 



(45) 



i=l 



where z^'^^ is the ith solution of z{t) in Eq. (44 1 



A. Quantum Control with Convoluted Pulses 

We start considering the case where e{t) is the re- 
sult of the convolution between a square pulse with a 
rise-time of 1.5 ns, a variable width, and a maximal am- 
plitude of 67.5 fieV, and a normalized Gaussian whose 
standard deviation is 3.7ns [see Fig. [l]^d)]. We plot 
in Fig. |6f^a)-(d) the corresponding singlet return prob- 
ability Ps as a function of iw and nominal final posi- 
tion Ef for different values of 71 and 72. We show in 
Fig. [6]Je)-(f ) traces taken at Sf = 3.985 meV (before the 
anti-crossing) and Sf = 3.95 meV. Here, we have cho- 
sen B = 100 mT, u = 4meV, r = 5/xeV and respec- 



tively 71 



10^ 



72 



10« 



[(a), blue trace (e)- 



(f)], 71 - lO^s-i, 72 = lO^s-i [(b), green trace (e)-(f)], 
71 = lQ^s-\ 72 = lO^s-i [(c), purple trace (e)-(f)], 
71 = lO^s-i, 72 = lO^s-i [(d), red trace (e)-(f)]. Our 
results indicate that charge noise strongly affects the dy- 
namics while spin relaxation only has a moderate effect. 
Although this behavior can be identified when comparing 
Fig.j6|^a) with Fig. [6](b) and Fig. IHJc), it is best seen in 
Fig.fmf) when comparing traces. An increase in the noise 
rate 72 leads to a substantial decrease of the oscillation 
visibility (blue and green traces). On the other hand, 
the visibility is only slightly diminished when relaxation 
is enhanced (blue and purple traces). 

Here, charge noise has two important consequences on 
spin dynamics. We have already mentioned that it leads 
to dephasing of the qubit when its state is a superposition 
of spin and charge states. The second consequence is a 
mechanism which competes with the LZSM interferom- 
etrjlSMse] gj^i^g ^^j^g hyperfine mediated anti-crossing is 
close to the charge anti-crossing for reasonable values of 
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Figure 6. (color online) Singlet return probability Ps as a 
function of the waiting time iw and final position et obtained 
for a Gaussian detuning pulse. Here, we have set B = 100 mT, 
u — 4meV, and r — 5/ieV for all plots. The parameters 71 
and 72 are (a) 71 = 10''s-\ 72 = 10** s'S (b) 71 = lO^'s'S 



72 



10^s-\ (c) 71 = 10^ s-^, 



72 



10*s-\ and (d) 71 



10^ s~^, 72 = 10'' s~^. Figures (e) and (f) show cuts for case 
(a) in blue, (b) in green, (c) in purple, and (d) in red. The 
cuts are respectively taken before the anti-crossing at et = 
3.985 meV and after at ei — 3.95 meV. These results show 
that charge noise is an important source of dephasing. 



B, charge noise also affects the dynamics during the pas- 
sage through the anti-crossing. As a result, the efficiency 
of the LZSM mechanism to create a coherent superposi- 
tion or to produce coherent interferences is hindered, and 
thus the optimal states populations are not reached. This 
behavior can be identified in Fig. when comparing 
traces obtained with same relaxation rate, but different 
charge noise rates (i.e. blue with green trace and purple 
with red trace). We clearly identify that the visibility of 
the oscillations is smaller for larger values of 72. 

It is not ideal to use detuning pulses with rise-times 
that are long compared to the dephasing time set by 
charge noise. Although we have chosen a trapezoid pulse 
with a rise-time of 1.5 ns, the effective rise-time due to 
the convolution is ss 10 ns. Under such conditions, the 
charge noise induced spin decoherence mechanism domi- 
nates for £ > £c- It thus becomes impossible to observe 
finite-time effects, as it can be seen in Fig.[6]Je). 

Energy relaxation processes have a weaker infiuence on 
spin dynamics than dephasing due to their dependence 
on the energy difference between the eigenstates of the 
system. These become only important when the system 
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is held in close vicinity of the anti-crossing, where the 
relaxation rate is at its maximum [c.f. Fig.pfa)]. 



B. Quantum Control with "Double Hat" Pulses 

We now consider a "double hat" detuning puls#^ [see 
Fig. [l|d)], whose leading and trailing edges are a se- 
quence of fast-slow-fast rise-time pulses as described in 
Sec. HU This control method has the benefit of achiev- 
ing a large adiabatic transition probability without ex- 
ponentially extending the overall pulse length. In terms 
of qubit manipulation, this means larger amplitude rota- 
tions can be achieved in one passage. 

We consider a detuning pulse s{t) with an amplitude 
of £f — Ei = 0.2 meV. The pulse reaches an amplitude of 
0.12 meV in 0.1 ns, it is then slowed down until it reaches 
an amplitude of 0.1325 meV, and finally it is brought 
to its maximal amplitude in 0.1ns. The rise time tsiow 
of the slow part of the pulse can be tuned freely. Here 
we compute Ps for tgiow = 6 ns. In Fig. [7]^a)-(d), we 
present the singlet return probability as a function of the 
final position ef and waiting time t^j for the same pre- 
viously considered parameters. The interference fringes 
are characterized, for e < ec, by three distinct regions 
showing an alternate oscillation amplitude for Ps- The 
darker regions coincide with detunings for which the pas- 
sage through the anti-crossing happened during one of 
the fast rise-times of e{t). Similarly, the bright region 
between e ~ 3.900 — 3.915 meV coincides with a passage 
through the anti-crossing with the slow rise-time portion 
of the detuning pulse. This region exhibits better oscil- 
lation visibility than what is obtained with a Gaussian 
pulse. These results are consistent with previous experi- 

This type of detuning pulse allows the direct obser- 
vation of finite-time effects in LZSM interferometry. We 
clearly identify coherent evolution of the qubit for e > Ec, 
Fig.[7[e) [blue and purple trace, 72 = 10^ s^^]. This indi- 
cates that it is possible to design complex pulses that can 
directly influence the competition between LZSM physics 
and charge noise. However, this is conditional on the de- 
phasing time scale associated with charge noise. We no- 
tice that if 72 = 10^ s~^, then it is impossible to identify 
any coherent behavior [green and red traces Fig. [7|[e)]. 
These results suggest that inhomogeneous dephasing due 
to nuclear spin fluctuations is not the only physical pro- 
cess that limits the coherence of a two-spin based qubit, 
but also charge noise plays a major role. This is an impor- 
tant result for future devices made out of Si/SiGe which 
are reaching a maturity level comparable to GaA 4^21581 
Silicon based devices are very interesting candidates for 
spin based quantum computing because the only stable 
isotope (^^Si) possessing a nuclear spin (/ — 1/2) has rel- 
atively low abundance, « 5%. Thus, hyperfine induced 
decoherence is weaker than in GaAs based nanostruc- 
ture^. 

Consequently, we distinguish three time scales that 
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Figure 7. (color online) Singlet return probability Ps as a 
function of the waiting time iw and final position ef obtained 
for a "double hat" detuning pulse. Pulse details are given in 
the text. Here, we have set B = 100 mT, u — 4meV, and 
T = 5 /ieV for all plots. The parameters 71 and 72 are (a) 
71 = 10'*s-\ 72 = W^s-\ (b) 71 = 10*s-\ 72 = W^s-\ 
(c) 71 = 10^s-\ 72 = lO^s-i, and (d) 71 = 10^s-\ 72 = 
10^ s~^. Figures (e) and (f) show cuts for case (a) in blue, (b) 
in dark green, (c) in purple, and (d) in dark red. The cuts are 
respectively taken before the anti-crossing at et = 3.985 meV 
and after at et = 3.908 meV. 



govern the physics of partial spin adiabatic passage in 
DQDs. There is the rise-time of the detuning pulse, the 
decoherence time T| , and associated with charge noise 
induced dephasing. Ideally, we would like to have rise- 
times shorter than to preserve any spin superposition 
state. But this would render adiabatic transitions un- 
likely and therefore seriously hinder manipulation of the 
qubit. Double hat pulses partially solve the problem, 
but for an even better result, it would be necessary to 
increase the coupling between the qubit states. In GaAs 
double quantum dots, it is possible to prepare a nuclear 
gradient field to enhance the hyperfine couplin^I^, with 
the advantage of extending T2*. In almost nuclear spin 
free systems (Si or G based DQDs), it is possible to use 
micro-magnets to artificially induce a coupling between 
the singlet state and triplet T+l^. 



C. Hadamard gate 

We now demonstrate that charge induced spin dephas- 
ing is not a fundamental roadblock in the implementation 
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Figure 8. (color online) It is possible to achieve Ps = 0.5 in 
the presence of charge dephasing by increasing the coupling 
between spin states to enhance the adiabatic transition prob- 
ability, (a) Singlet state return probability for a prepared nu- 
clear state with tsiow ~ 4ns, 71 = 10^ s~^, and 72 — 10* s~^. 
(b) In case of fastest charge decoherence, it is still possible to 
achieve perfect beam splitting by further increasing the cou- 
pling. This is achieve by increasing tsiov = 8 ns. We have 
71 = 10*s-\ and 72 = 10^ 8"^ 



of high fidelity quantum gates using LZSM interferome- 
try. We specifically show that it is possible to achieve 
Ps = 0.5, this particular value being a signature of a 
perfect 50 : 50 beam splitter. For quantum information 
processing, this allows a Hadamard gate Eq. Q, w hich 
is required to perform some quantum algorithrna^Zl. 

In the following, we consider a nuclear state that has 
been prepared according to the method developed in 
Ref. [m i.e. a nuclear gradient field has been generated. 
We model such a prepared nuclear state by considering 
a mean difference field |(5i3n| |-B„,2 — -Bn,i| ~ lOmT. 
The Gaussian distribution of nuclear fields Eq. (25) is 
changed to reproduce the mean SB^. To take into ac- 
count reduced nuclear fluctuations we consider a stan- 
dard deviation tTprep = cr/6. We present in Fig.|8[a) a cut 



of the return singlet probability for a "double hat" pulse 
along £f = 3.911 meV with ts\ow = 4ns, B = lOOmT, 
u = 4meV, t — 5^eV, 71 = 10^ s~^, and 72 — 10^ s^^. 
The results demonstrate that for this particular choice of 
parameters, it is possible to achieve Ps = 0.5. 

The remaining question is whether it is still possible to 
achieve Ps — 0.5, if the charge noise rate 72 is increased. 
As expected from our previous analysis, and for the pre- 
vious set of parameters, it is not possible. However, 
if we increase the coupling by increasing t^iow, we can 
again achieve perfect beam splitting. This is presented 



in Fig.[8](b), where we have plotted a cut of the singlet re- 
turn probability along e{ = 3.9075 meV with igiow = 8 ns, 
B = lOOmT, u = 4meV, r = bfieV, 71 = 10^ s~^, and 
72 = 10^ s^^. Although these results demonstrate that in 
principle it is possible to reach Ps = 0.5, it is still neces- 
sary to study the average gate fidelity. This goes beyond 
the scope of this paper and will be presented elsewhere. 



V. CONCLUSIONS 

We have developed a master equation formalism to 
study the dynamics of a S-T_|_ qubit in the vicinity of 
the hyperfine mediated anti-crossing. In comparison with 
previous theories that only included decoherence due to 
nuclear spin fluctuation, we also include spin relaxation 
(e.g. due to spin-orbit) and spin dephasing due to charge 
dynamics phenomenologically. We have also derived a 
new effective two-level spin-charge Hamiltonian. In addi- 
tion to previous theories, we take into account the charge 
degree of freedom. This property originates from the 
ground state singlet |S) being a superposition of different 
charge configurations. Although the effective Hamilto- 
nian is strictly not of the LZSM form, the peculiar form 
of the coupling constant allows to still reason in terms of 
LZSM physics. 

Using our formalism, we have calculated the singlet 
return probability for two types of detuning pulses. For 
each case, we have compared results for different values 
of 71 and 72. Our findings suggest that LZSM spin in- 
terferometry can be inhibited by charge dynamics, and 
thus charge coherence has to be treated on equal foot- 
ing with spin coherence. We have indeed demonstrated 
that the visibility of the spin qubit coherent oscillations 
is sensitive to the time scale associated with charge in- 
duced spin dephasing. We have also shown that this in- 
terplay between charge and spin can prevent the observa- 
tion of finite-time oscillations. Interestingly, this last fact 
could lead to the development of a measurement protocol 
that would allow to quantify the time scale associated to 
charge decoherence. 

Finally, we have demonstrated that in the case of a 
prepared nuclear state with a mean gradient field and a 
"double hat" detuning it is possible to achieve Ps = 0.5. 
This demonstrates that an equally weighted superposi- 
tion of S and T_|_ can be reached. This is an important 
result for quantum information processings as it corre- 
sponds to the achievement of a Hadamard gate. It is 
an essential element in the implementation of important 
quantum algorithms like Shor's algorithm for period find- 
ing, and which ultimately allows to factorize large inte- 
gers. Although we have supposed that a gradient field is 
created via the method developed in Ref.[T31 it is possible 
to prepare the nuclear state via weak measurements as 
proposed in Ref. |49] to freeze the nuclear spin dynamics 
and create an artificial splitting by means of a micro- 
magnell^. 
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